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Abstract 

Using a relationship between the moments of the probability dis- 
tribution of times between the two consecutive trades (intertrade time 
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£^ . distribution) and the moments of the distribution of a daily number of 

trades, we show that the underlying point process is essentially non- 
Markovian. A detailed analysis of all trades in the EESR stock on 
the Moscow International Currency Exchange in the period January 
2003 - September 2003, including correlation between intertrade time 
CN ' intervals is presented. A power-law decay of the correlation function 

provides an additional evidence of the long-memory nature of the se- 
, ries of times of trades. A data set including all trades in Siemens, 

Commerzbank and Karstadt stocks traded on the Xetra electronic 
^ ■ stock exchange of Deutsche Boerse in October 2002 is also considered. 
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The new field of econophysics, in particular part of these activities related 
to the description of financial markets, has drawn considerable attention in 
recent years jUEUHj. One of the central issues under discussion is understand- 
ing fundamental origins of the stock price formation process. The progress 
in this direction is very rapid, see, e.g., [3J El El HI El El ED! and references 
therein. The ultimate goal of the quantitative approach to the description of 
stock price dynamics is constructing a microscopic theory, in which this dy- 
namics originates from the agents' activity in the financial market generating 
patterns of demand - supply interaction, resulting in price formation. At a 
more phenomenological level the description of stock price evolution focuses 
on the dynamics of accomplished trades. 

Among the most important characteristics of financial markets is their 
activity. On the trade-by-trade basis one can think of two possible ways of 
characterizing it and thus setting the clock measuring the operational time. 
The first possibility is to follow the volume traded fTT| . the second is to an- 
alyze the time-dependent frequency of trading operations. Below we shall 
follow the latter route. The importance of effects related to the varying trad- 
ing frequency is well recognized. In particular, in 12J it was suggested, that 
they significantly influence the expected return during intensive speculative 
trading. In [4 j it was also argued that the empirically observed long-range 
correlations in trading frequency [T31 H] directly induce the observed long- 
range correlations of volatility. Let us also mention a related topic of great 
practical interest - a study of seasonality effects (e.g., intraday and intraweek 
activity patterns) and of the corresponding tuning of trading strategies, see 
e.g. [H. 

Below we shall study the trade dynamics at the smallest level of resolu- 
tion, focusing ourselves on the properties of a (stochastic) process generating 
the times of trades. Given a description of such directing point process T(t) 
generating the times of trades ti < £2 , the temporal evolution of quanti- 
ties of interest such as price or volume is that of a subordinated stochastic 
process JE] P{T{t)). The form the corresponding evolution equation for the 
probability distributions of, e.g., price fl(P,t) and its characteristic features 
are thus decisively dependent on the properties of the directing point process 
T(t). Generically one can distinguish the following possibilities: 

• Type 1. The time at which the n-th point (trade) takes place is 
completely independent of the times at which the previous n — 1 points 
were generated. In this case one deals with the Poisson distribution 
of the number of trades in any fixed time interval. This point process 
is fully characterized by an exponential probability distribution for the 
intertrade time r — t n — t n _i, ip(r) = 7~o~ exp(— t/tq). The evolution of 
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the probability density Q(P,t) is described by a differential (in time) 
equation. 

• Type 2. The time of the n-th trade is correlated with that of the 
previous, n — 1-th, trade, so that the the point process has a unit mem- 
ory depth (and is thus still Markovian). The point process T(t) is still 
fully characterized by some non- exponential intertrade time probability 
distribution i[>(t), but the evolution equation for Q(P,t) is no longer 
differential, but an integral one, of continuous-time random walk type 
[To] 3 . The correlation between the number of trades in non-overlapping 
time windows can in principle be computed [T7], but the resulting ex- 
pressions are quite involved. 

• Type 3. The time of the n-th trade depends on r > 1 times of the 
previous trades t n _i, t n _ 2 , t n _ r . In this case the point process is a 
long-range memory non-Markovian one with a memory depth equal to 
r, and the corresponding evolution equation for Q(P, t) is a complicated 
integral one with a kernel depending on r — 1 arguments. If all temporal 
scales are effectively involved, r is infinite. 

The main question is, of course, whether the times of trades correspond to 
a long-range memory process with past events strongly influencing the new 
ones. Let us first consider a standard, although not very reliable, measure 
of the amount of memory present in the series - the (normalized) correlation 
function of the intertrade time intervals 



The correlation is computed by first averaging over all intervals separated by 
a fixed number of intervals d for a given day and then averaging over all days. 
The correlation (0) for all trades, from January 2003 till September 2003, in 
the EESR stock traded on the MICEX stock exchange is shown, together 
with the errorbars characterizing its variation on the day-to-day basis, in 
Fig. ^ The correlation function C T (d) shown in Fig. Q]is clearly powerlike, 
C T (d) ~ 0.15/d ' 29 , thus demonstrating a behavior typical for long-range 
processes. 

To get a somewhat more precise picture of the memory - related pattern 
of the set of times of trades, let us consider a relation between the moments 
of the distribution pat(AT) of a number of trades in some given time interval 

3 The application of CTRW to the analysis of financial data was discussed in a number 
of papers QB HJ H3 EU HI E3 El 



C T {d) 





3 



EESR 200301—200309 




Distance in ticks 



Figure 1: Normalized correlation of intertrade time intervals in tick time. 
Solid line shows a powerlike fit of the form 0.15/<i ' 29 . Upper and lower lines 
show the standard deviation of the measured correlation characterizing its 
variation on the day-to-day basis. 



AT and those of the intertrade time distribution ip(r) that holds for the unit 
memory depth point process in the large - AT limit ^Zj. More specifically, 
for the processes of Type 1 and Type 2 there exist, for a sufficiently large time 
interval AT, a relation between the first two moments of the distribution of 
the number of trades within this interval (./Vat) and (N AT ) and the first two 
moments of the intertrade time distribution = (r) and ip2 = (t 2 ) 

(N AT )-(N AT ) — —j— = (N A t) (2) 

Corrections to the relation (J2J) are of order 1/AT. Equivalently, 

„ _ (Nj T ) - (Nat) 2 ^ zA - n (*\ 

so that if the underlying point process T(t) is indeed of Type 1 or Type 2, 
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one should have Pn/Pt — 1- Let us note, that in the simplest poissonian 
case pn = p T = 1. Let us also stress that provided the underlying point pro- 
cess has unit memory length, the equation (PJ) holds for any intertrade time 
probability distribution ip{r) having a finite second moment ip2- Derivation 
of Eq. (JHJ) is sketched in the Appendix. 

Equation Q provides a straightforward possibility of establishing the 
nature of the point process T(t) by computing the quantities in its left- 
and right-hand side from an observed distribution of a number of trades in 
some fixed interval AT and a corresponding underlying distribution of the 
intertrade times. In particular in the case of the ensemble of independent 
points (Poisson case) p^ = p T = 1, in the case of the CTRW - type process 
having unit memory depth a more general relation of Eq. © holds and, 
finally, for the long-range memory process the relation in Eq. (jSJ) should 
break down. 

In the data set used for our analysis the interval AT corresponds to 
one trading day (10.30 - 18.45 on the MICEX exchange and 9.00 - 20.00 
on the Xetra electronic exchange of Deutsche Boerse in Frankfurt), and the 
considered ensemble of trades consists of all trades in EESR in January - 
September 2003 (MICEX) and Siemens, Commerzbank and Karstadt (liquid, 
medium liquid and relatively illiquid stock) in October 2002 on Xetra. The 
results are summarized in Table 1: 



Table 1 



Und 


(N AT ) 


(r) (seq) 


Pr 


Pn 


Pn/Pt 


EESR 


6823 


4.5 


10.3 


492.8 


109.4 


SIE 


4364 


9.1 


3.8 


124.8 


32.8 


CBK 


1856 


21.3 


3.1 


274.4 


88.5 


KAR 


373 


104.1 


4.5 


39.7 


8.8 



From the results presented in Table 1 it is clear, that poissonian Type 
1 and CTRW Type 2 processes are excluded as candidate point processes 
generating the times of the trades. This leaves us with the only remaining 
possibility of Type 3 non-Markovian long-range memory process. This con- 
clusion is in obvious agreement with the long-range correlation between the 
intertrade time intervals for EESR shown in Fig. (JIJ. 

To provide some statistical backing to this conclusion, we have generated 
the times of trades using the empirically observed distribution over intertrade 
times as the EESR data. The thus generated set of points is characterized 
by the observed distribution, but is otherwise uncorrelated 4 . The considered 

4 The distributions over intertrade time intervals were studied in a number of papers 
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number of trade intervals AT was equal to the number of trading days in the 
EESR data. This simulation effectively reconstructs a unit memory depth 
CTRW process. The results of this simulation are shown in Table 2: 



Table 2 



Und 


(N AT ) 


(r) (seq) 


Pr 


Pn 


Pn/Pt 


EESR (data) 


6823 


4.5 


10.3 


492.8 


109.4 


EESR (sim.) 


6816 


4.5 


12.4 


8.2 


0.66 



We see, that although due to statistical limitations (finite size corrections, 
etc.) the simulation does not reproduce the theoretical value of Pn/Pt = 1, 
it is sufficiently close to it, while at the same time missing the experimental 
ratio by two orders of magnitude. 

Conclusion 

We conclude that the point process generating the times of trades is of 
long-range memory non-Mar kovian nature. This excludes random walk and 
continuous-time random walk (CTRW) processes as models describing high 
frequency stock price dynamics in real time. 
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Appendix A Derivation of Eq. (3) 

In this Appendix we sketch, following [T7j, the derivation of Eq. (J3J). The 
probability of having A" trades in the interval [0,T] is conviniently written 
as a convolution of a probability of having A" trades in the interval [0, t] (the 
AT-th trade falling into the infinitesimal vicinity of t) w^ft) and a probability 
of having no trades in the remaining interval [t, T] wq{T — t)\ 

p N {T) = ( T dtw (T-t)w N (t) (A.l) 
Jo 

so that the Laplace transform of Pn(T) reads 

p N (X) = J dTe~ XT p N (T) = w (\)w N (\) = 1 ~^ X) j>(\) N - 1 , (A.2) 

where we took into account that in the considered unit memory depth case 
all probabilities in question are fully determined by the intertrade (waiting) 
time distribution tp(r), so that 5 w)at(A) = w(A) Ar_1 and, for w we have, 
recalling that Wo (r) = dr' 4>(t'), Wq(\) = (1 — ^(A))/A. 

To compute the moments of trade's multiplicity distribution it is conve- 
nient to introduce a generating function 

oo oo nk 

G{9\T) = £e fc V0O = 1 + E -u(N k ) (A.3) 
k=i k=i K - 

5 Here we assume that there is always a first trade at the beginning of the time interval 
in question 
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so that the moments are given by its derivatives at the origin 

(-) 

Combining Eqs (jA.2|A.3() one gets the following expression for the gener- 
ating function: 

Equation (j^J follows now by substituting the expansion of if>(\) in the mo- 
ments of ip{r) 

V»(A) = / rfre- A ^(r) = 1 - ^ A + - ^ 2 A 2 + ... (A.6) 
Jo 2 



into Eq. (jA.5|) and computing, through Eq. (jA.4j) . the leading contributions 
to (N T ) and (iV|) in the large T limit 
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